APPARATUS FOR AND METHOD OF MEASURING CLOCK SKEW 



BACKGROUND OF THE INVENTION 

The present invention relates to a clock skew measurement apparatus and 
5 a clock skew measurement method by which a skew between a plurality of clock 
signals distributed by a clock distribution circuit is measured. 

A clock skew has conventionally been estimated statistically using a time 
interval analyzer or a frequency counter. That is, as shown in Fig. 1, for 
example, when a reference clock signal CLK g from a clock signal source 11 is 

10 distributed and supplied to registers 12j and 12k, each of the clocks under 

measurement CLKj and CLKk in the respective registers 12j and 12k is inputted 
to a time interval analyzer 13, where a timing difference between a zero-crossing 
point of the clock signal under measurement CLKj and a zero-crossing point of 
the clock signal under measurement CLKk is measured, and a fluctuation of the 

15 timing difference is measured by a histogram analysis. An example of clock 
skew measurement using a time interval analyzer 13 is described in, for example, 
"Jitter Analysis Clock Solutions", Wavecrest Corp., 1998. 

However, there is a problem in this clock skew measurement method 
using the time interval analyzer that it takes a long time to acquire a number of 

20 data required for the histogram analysis since, it requires an intermediate 

dead-time until clock skew measurement can be repeated. In addition, in the 
clock skew measurement method using the time interval analyzer, a skew 
between clocks having different frequencies cannot be measured. Therefore, a 
new clock skew measurement method is required for a precise control of a local 

25 clock and a global clock. 
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It is an object of the present invention to provide a measurement 
apparatus of a skew between clocks and its method by which a skew between 
clocks can be measured in a short time compared with a conventional apparatus 
and its method. 

5 It is another object of the present invention to provide a measurement 

apparatus of a skew between clocks and its method by which a skew between 
clocks having different frequencies can be estimated. 

SUMMARY OF THE INVENTION 

10 The clock skew measurement apparatus according the present invention 

comprises timing jitter estimation means for estimating timing jitter sequences of 
a plurality of clock signals under measurement, and clock skew estimation 
means to which the plurality of timing jitter sequences are inputted for 
calculating a timing difference sequence between the plurality of timing jitter 

15 sequences to output a clock skew sequence. 

In addition, it is desirable that the clock skew measurement apparatus 
includes second clock skew estimation means to which the clock skew sequences 
are inputted for obtaining a difference between the plurality of clock skew 
sequences. 

20 In addition, it is desirable that the clock skew measurement apparatus 

includes frequency multiplication means to which the timing jitter sequence is 
inputted for multiplying a frequency of the timing jitter sequence to output a 
timing jitter sequence of the frequency-multiplied clock signal under 
measurement. 

25 In addition, it is desirable that the clock skew measurement apparatus 
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includes deterministic clock skew estimation means for estimating a timing error 
between ideal clock edges of the plurality of clock signals under measurement to 
output a deterministic component of clock skew. 

In addition, it is desirable that the clock skew measurement apparatus 

5 includes clock skew detection means to which the clock skew sequence is 
inputted for obtaining clock skew values between the clock signals under 
measurement from the clock skew sequence. 

In addition, it is desirable in the clock skew measurement apparatus that 
the timing jitter estimation means comprises analytic signal transformation 

10 means for transforming a clock signal under measurement into a complex 
analytic signal, instantaneous phase estimation means for obtaining an 
instantaneous phase of the analytic signal, linear phase removing means for 
removing a linear instantaneous phase from the instantaneous phase to obtain an 
instantaneous phase noise, and zero-crossing sampling means to which the 

1 5 instantaneous phase noise is inputted for sampling only the instantaneous phase 
noise data at timings close to zero-crossing timings of a real part of the analytic 
signal to output a timing jitter sequence. 

In addition, it is desirable in the clock skew measurement apparatus that 
the analytic signal transformation means comprises bandpass processing means 

20 to which the clock signal under measurement is supplied for extracting only 
components around a fundamental frequency from the clock signal under 
measurement to band-limit the clock signal under measurement, and Hilbert 
transformation means for Hilbert-transforming an output signal of the bandpass 
processing means to generate a Hilbert pair of the input signal. 

25 In addition, it is desirable in the clock skew measurement apparatus that 
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the analytic signal transformation means comprises time domain to frequency 
domain transformation means to which the clock signal under measurement is 
supplied for transforming the clock signal under measurement into a both-sided 
spectrum signal in frequency domain, band-limiting means for extracting only 
5 components around a positive fundamental frequency of the clock signal under 
measurement in the both-sided spectrum signal, frequency domain to time 
domain transformation means for inverse-transforming an output of the 
band-limiting means into a signal in time domain, 
y In addition, it is desirable in the clock skew measurement apparatus that 

pf 10 the analytic signal transformation means comprises a buffer memory to which 
^ the clock signal under measurement is supplied for storing therein the clock 
pU signal under measurement, means for extracting the signal in the sequential order 
s n from the buffer memory such that the signal being extracted is partially 
Si overlapped with the signal extracted just before, means for multiplying each 
W 15 extracted partial signal by a window function, means for transforming each 
1-* partial signal multiplied by the window function into a both-sided spectrum 

signal in frequency domain, band-limiting means for extracting only components 
around a positive fundamental frequency of the clock signal under measurement 
from the transformed both-sided spectrum signal in frequency domain, means for 
20 inverse-transforming an output of the band-limiting means into a signal in time 
domain, and means for multiplying the transformed signal in time domain by an 
inverse number of the window function to obtain a band-limited analytic signal. 

In addition, it is desirable that the clock skew measurement apparatus 
includes AD conversion means to which the clock signal under measurement is 
25 inputted for digitizing an analog signal to convert it to a digital signal. 
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In addition, it is desirable that the clock skew measurement apparatus 
includes waveform clipping means to which the clock signal under measurement 
is inputted for removing amplitude modulation components of the clock signal 
under measurement to extract only phase modulation components of the clock 
5 signal under measurement. 

In addition, it is desirable in the clock skew measurement apparatus that 
the clock skew detection means is peak-to-peak detection means for obtaining a 
difference between the maximum value and the minimum value of the supplied 
clock skew sequence 

10 In addition, it is desirable in the clock skew measurement apparatus that 

the clock skew detection means is RMS detection means for obtaining a 
root-mean-square value (RMS value) of the supplied clock skew sequence. 

In addition, it is desirable in the clock skew measurement apparatus that 
the clock skew detection means is histogram estimation means for obtaining a 

1 5 histogram of the supplied clock skew sequence. 

In addition, it is desirable in the clock skew measurement apparatus that 
a pass bandwidth of the clock signal under measurement is variable in the 
analytic signal transformation means. 

In addition, it is desirable in the clock skew measurement apparatus that 

20 the timing jitter estimation means further includes low frequency phase noise 
removing means to which the instantaneous phase noise is inputted for removing 
low frequency components of the instantaneous phase noise to output the 
instantaneous phase noise from which the low frequency components have been 
removed to the zero-crossing sampling means. 

25 In addition, the clock skew measurement method according to the present 
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invention includes a step of estimating timing jitter sequences of a plurality of 
clock signals under measurement, and a step of calculating a timing difference 
sequence between the plurality of timing jitter sequences to estimate a clock 
skew sequence. 

5 In addition, it is desirable that the clock skew measurement method 

includes a step of obtaining a difference between the plurality of clock skew 
sequences, to estimate a clock skew sequence. 

In addition, it is desirable that the clock skew measurement method 
includes a step of assigning the timing jitter values of the timing jitter sequence 

10 to provide a timing jitter sequence that corresponds to the frequency-multiplied 
clock signal under measurement. 

In addition, it is desirable that the clock skew measurement method 
includes a step of estimating a timing error between ideal clock edges of the 
plurality of clock signals under measurement to estimate a deterministic 

1 5 component of clock skew. 

In addition, it is desirable that the clock skew measurement method 
includes a step of obtaining clock skew values of the clock signals under 
measurement from the clock skew sequence. 

In addition, it is desirable in the clock skew measurement method that the 

20 step of estimating a timing jitter sequence includes a step of transforming a clock 
signal under measurement into a complex analytic signal, a step of obtaining an 
instantaneous phase of the clock signal under measurement from the analytic 
signal, a step of removing a linear instantaneous phase from the instantaneous 
phase to obtain an instantaneous phase noise, and a step of sampling, using the 

25 instantaneous phase noise as an input, only the instantaneous phase noise data at 
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timings close to zero-crossing timings of a real part of the analytic signal to 
output a timing jitter sequence. 

In addition, it is desirable in the clock skew measurement method that the 
step of transforming the clock signal under measurement into an analytic signal 
5 includes a step of extracting only components around a fundamental frequency 
from the clock signal under measurement to band-limit the clock signal under 
measurement, and a step of Hilbert- transforming an output signal of the 
band-limiting processing means to generate a Hilbert pair of the input signal. 

In addition, it is desirable in the clock skew measurement method that the 

10 step of transforming the clock signal under measurement into an analytic signal 
includes a step of transforming the clock signal under measurement into a 
both-sided spectrum signal in frequency domain, a step of extracting ordy 
components around a positive fundamental frequency of the clock signal under 
measurement in the both-sided spectrum signal, and a step of 

15 inverse-transforming an output of the band-limiting means into a signal in time 
domain. 

In addition, it is desirable in the clock skew measurement method that the 
step of transforming the clock signal under measurement into an analytic signal 
includes a step of storing the clock signal under measurement in a buffer 

20 memory, a step of extracting the signal in the sequential order from the buffer 
memory such that the signal being extracted is partially overlapped with the 
signal extracted just before, a step of multiplying each extracted partial signal by 
a window function, a step of transforming each partial signal multiplied by the 
window function into a both-sided spectrum signal in frequency domain, a step 

25 of extracting only components around a positive fundamental frequency of the 
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clock signal under measurement from the transformed both-sided spectrum 
signal in frequency domain, a step of inverse-transforming the band-limited 
spectrum signal into a signal in time domain, and a step of multiplying the 
transformed signal in time domain by an inverse number of the window function 
5 to obtain a band-limited analytic signal. 

In addition, it is desirable in the clock skew measurement method that the 
step of estimating a deterministic component of clock skew between the clock 
signals under measurement includes a step of using the instantaneous linear 
phases of the plurality of clock signals under measurement as inputs and 

10 obtaining a difference between initial phase angles of the instantaneous linear 
phases to obtain a deterministic component of clock skew. 

In addition, it is desirable in the clock skew measurement method that the 
step of estimating a deterministic component of clock skew between the clock 
signals under measurement includes a step of using the timing jitter sequences of 

15 the plurality of clock signals under measurement as inputs and obtaining a 
correlation between the timing jitter sequences to estimate strongly correlated 
clock edges and to obtain an offset value that relates two clock edge with each 
other. 

In addition, it is desirable in the clock skew measurement method that the 
20 step of estimating a deterministic component of clock skew between the clock 
signals under measurement includes a step of using the plurality of clock signals 
under measurement as inputs for obtaining a mean value of zero-crossing timing 
errors between the clock signals under measurement to obtain a deterministic 
component of clock skew. 
25 In addition, it is desirable that the clock skew measurement method 
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includes a step of performing a waveform clipping of the clock signal under 
measurement in order to remove amplitude modulation components of the clock 
signal under measurement and extract only phase modulation components of the 
clock signal under measurement. 
5 In addition, it is desirable in the clock skew measurement method that in 

the step of obtaining the clock skew, a difference between the maximum value 
and the minimum value of the clock skew sequence is obtained to calculate a 
peak-to-peak value. 

In addition, it is desirable in the clock skew measurement method that in 
10 the step of obtaining the clock skew, a root-mean-square value of the clock skew 
sequence is obtained to calculate an RMS value. 

In addition, it is desirable in the clock skew measurement method that in 
the step of obtaining the clock skew, a histogram data of the clock skew 
sequence is obtained. 

15 In addition, it is desirable in the clock skew measurement method that the 

step of estimating the timing jitter further includes a step of using the 
instantaneous phase noise as an input for removing low frequency components 
of the instantaneous phase noise. 

The principle of the present invention will be explained. In this 

20 explanation, a clock signal in a microprocessor unit is used as a clock signal 
under measurement. 

CLOCK SKEW MEASUREMENT METHOD 

First, a clock skew will be defined. As shown in Fig. 2, a clock skew is 
25 given by the difference between delay times xj and T cd k of the clock signals 
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CLKj and CLKk to arrive at the respective registers 12j and 12k when the clock 
signals CLKj and CLKk are derived from a reference clock signal CLK g of a 
clock signal source 11 as a reference time point of, for example, a clock 
distribution network. That is, the difference is given by an equation (1). 

5 

4"w( nT ) = ^d( nT )-^( nT ) (D 

Fig. 3 shows by dotted lines a reference clock signal CLK g having a 
fundamental period T, and clock signals CLKj and CLKk. Also shown in Fig. 3 
10 are a difference t c j (nT) between a rising edge of the reference clock CLK g and a 
rising edge of the clock signal CLKj, a difference icd k (nT) between the rising 
edge of the reference clock CLKg and a rising edge of the clock signal CLKk 
(n=0, 1, 2, ••■), and a difference between the icd" (nT) and i C d k (nT), i.e., a clock 
skew T ske w' k (nT). 

15 As shown in Fig. 4, it is assumed that rising edges of the clock signals 

CLKg, CLKj and CLK k occur at time t C d 8 (nT), t^ (nT) and t C d k (nT) respectively, 
and ideal clock edges (clock edges correspond to jitter-free clock signals) of the 
clock signals CLK g , CLKj and CLKk rise at time (nT) g , (nT)j and (nT)k 
respectively. Then the delay time z c j (nT) or i C d k (nT) between the clock 

20 departure of the clock signal CLKj or CLKk and the arrival at the respective 
registers 12j or 12k is expressed by following equation. 
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^id ( nT > = tid C^T) - tf d (nT) 

= ^^-(nT^J-^Cn^-TOgl+Kn^j -(nT) g } 



x Skew + 



A<|) J [n] 



-A()) g [n] 



\2%j 



[sec] 



(2) 



' c cd( nT ) = t cd( nT )- t ^( nT ) 

= [t^ d (nT)-(nT) k ]-[t| i (nT)-(nT) g ]+{(nT) k -(nT) g } 



- T 8 ' k + 

" T Skew + 



A(t> k [n] 



\2%j 



-A(|) g [n] 



v 27t yj 



[sec] 



(3) 



5 In this case, 



10 



15 



4k J ew=( nT )j-( nT )g 



and 



ikew=( nT )k-( nT )g 



[sec] 



[sec] 



(4) 



(5) 



are a time difference between an ideal clock edge of the clock signal CLKj and 
an ideal clock edge of the reference clock signal CLK g5 and a time difference 
between an ideal clock edge of the clock signal CLK k and an ideal clock edge of 
the reference clock signal CLK g , respectively, and those time differences 
correspond to deterministic components of clock skew (deterministic clock skew 
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values) that are determined by their respective paths. In addition, A<J> g [n](T g /27t) 
Hd § (nT)-(nT) g ), Atf[n](y2%) (=*J (nT)-(nT)j) and A(j) k [n] (Tk/27i)(=tcd k 
(nT)-(nT)k) express timing jitter sequences (each unit is second) of the clock 
signals CLK g , CLKj and CLKt, respectively. When the equations (2) and (3) 
are substituted in the equation (1), a clock skew T S kew between the clock signal 
CLKj and the clock signal CLKk is expressed by the following equation. 



T 8 ' k + 

T Skew + 



A<|> k [n] 



-A<j) 8 [n] 



X 



T 8 ' J + 

x Skew + 



A<H[n] -i _ A ((,g[n] -5- 

^ZTT J \ J 



[sec] (6) 



_T Skew + 



A^ k [n/il-A(|)j[n/i- N 

\Z% J \^Z% j 



10 The second term 

f^v \ . fT-Y 
A* k [n] -A(H[n] -i- 
\Zit j yz% J 

of the equation (6) corresponds to a random spread (random components) of 
clock skew caused by timing jitter of each clock signal. Then in the clock skew 
measurement method according to the present invention, an random distribution 
15 is obtained by combining quantities each representing how much a clock edge of 
each clock signal is deviated from an ideal clock edge, i.e., timing jitters of the 
clock signals. In this case, it is assumed that the fundamental period of the 
distributed clock signal CLKj is generally equal to the fundamental period of the 
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distributed clock signal CLKk (Tj=Tk). In addition, 

4kew=( nT )k-( nT )j M ( ? ) 



10 



is a difference between a rising edge of an ideal clock of the clock signal CLKj 
and a rising edge of an ideal clock of the clock signal CLKk, and this Tskew ' k is a 
deterministic component of clock skew determined by the paths of clock 
distribution network. 

A deterministic clock skew value TskeJ* can be obtained by, for example, 
obtaining instantaneous phases of corresponding two clock signals under 
measurement CLKj and CLK k , and then by calculating a difference between their 
linear phase components. Fundamental components of the clock signals CLKj 
and CLKk are expressed as cosine waves by following equations (8) and (9), 
respectively. 



15 Xj(t)=AjCos((|)j(t)) = AjCos 



T 



t + (|> J 0 -A(|>J(t) 



(8) 



Xfc (t) = Ak cos((|> (t)) = Afc cos 



W T k; 



t + ^-A^Ct) 



(9) 



Here, each of instantaneous phases (^(t) and §\t) of the respective Xj(t) 
20 and Xk(t) is expressed by a sum of a linear instantaneous phase component 27rt/TL 
(L=j, k) corresponding to a fundamental period T L (L=j, k), an initial phase angle 
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(|>o (L=j, k), and an instantaneous phase noise component A(|) L (t) (L=j, k). 



* J (t) = 



1 ' 2%^ 
T 



t + (|) J 0 -A<(»j(t) 



[rad] 



(10) 



r(t)= 



f 2k^ 



t+^-A(|) k (t) 



[rad] 



(11) 
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An estimation method of an instantaneous phase of a clock signal will be 
explained later. In the equations (10) and (11), when A<()(t)=0 is assumed, linear 
instantaneous phases of the respective jitter-free clock signals CLKj and CLKk 
can be estimated by following equations. 



f \ 
2% 



T 



[rad] 



(12) 



15 



20 



^linear ( t ) = 



t+«t>o 



[rad] 



(13) 



In this case, the ideal rising edge of the respective clock signals CLKj and CLKk 
occur at time (nT)j or (nT)k when each linear instantaneous phase of the left side 
becomes (2n7c-7t/2). And the following relations can be obtained from the 
equations (12) and (13), respectively. 
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(nT) i = 



2m- 



\2kJ 



[sec] 



(14) 



(nT) k = 



2n7i- 



(3-*. 



27T 



[sec] 



(15) 



5 Therefore, the equations (14) and (15) are substituted in the equation (7) to 
obtain the deterministic clock skew value iskew j,k by the following equation. 



4kew=( nT )k-(nT)j 



= 2n%-\ — 



2% 



2mt-f ^ 



-* J o 



\2%j 



= (<!>o-<t>o) 



kxTo 



2tt 



[sec] (16) 



10 In general, the fundamental periods of the respective clock signals CLKj and 
CLKk are the same (Tj=Tk). That is, a deterministic clock skew value between 
the two clock signals under measurement can be obtained as a difference 
between initial phase angles <|>o and <|>o k in the linear instantaneous phases of the 
respective two clock signals under measurement. 

15 Here, an initial phase angle <|>o of a clock signal under measurement can 

be obtained by applying a linear line fitting by least squares method to an 
instantaneous phase <|)(k) ( k is a digitizing time point), and by selecting $ 0 so 
that a formula (17) becomes minimum. 
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2((t>(k)-(«ok+^o)) 2 (17) 
k=l 

In this case, an initial phase angle to be obtained is given by the following 
equation (18). 



2N(2N + l)Z<Kk)-6£k(Kk) 

♦o = — — ( 18 ) 

Vo N(N-l) v } 



A parameter estimation by a linear line fitting is described in, for example, 
"Random Data: Analysis and Measurement Procedure" by J.S. Bendat and A.G. 

10 Piersol, 2nd ed., John Wiley & Sons, Inc., p.362, 1986. 

In addition, an initial phase angle <|>o of a clock signal under measurement 
x(t) can be obtained by applying a cosine wave fitting by least squares mthod to 
a clock waveform data x(k) or its fundamental sine wave component, and by 
estimating ^ 0 using a maximum likelihood estimation method so that a formula 

15 (19) becomes minimum. 



N f 

^ x(k)-Acos 



k+< 



AY 



J) 



(19) 
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In this case, an initial phase angle to be obtained is given by the following 
equation (20). 
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= -arctan 



^ N 2% ^ 

£x(k)sin— k 

k=l i_ 

N 2tc 
5>(k)cos— k 

Vk=i 1 



(20) 
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A parameter estimation by a maximum likelihood estimator is described in, for 
example, "Fundamentals of Statistical Signal Processing: Estimation Theory" by 
S.M. Kay, Prentice-Hall Inc., pp.167-172, 1993. 

In the above description, it is assumed that corresponding clock edges of 
the two clock signals under measurement are not spaced apart from each other 
equal to or more than one period. When the corresponding clock edges of the 
two clock signals under measurement are spaced apart from each other equal to 
or more than one period, a deterministic clock skew value is given by a sum of 
the difference between the initial phase angles and an offset time noffsetTo of a 
clock edge, i.e., by the following equation (21). 



15 



20 



T ^ = (*i"*o)^ +n offsetT 0 



[sec] 



(21) 



A clock signal distributed from the clock signal source has a strong 
cause-and-result relationship with a clock signal of the clock signal source. As 
a result, in general, a phase noise (timing jitter sequence) of a distributed clock 
signal shows a tendency similar to that of a phase noise (timing jitter sequence) 
of the clock signal source. For this reason, timing jitter sequences A^n] and 
A(|) k [n] of the respective clock signals CLKj and CLK k distributed from the same 
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clock signal source show mutually similar tendencies as shown in, for example, 
Fig. 5 A and Fig. 5B. Therefore, an amount of offset time n 0 ff S et of the 
corresponding clock edges of the respective clock signals CLKj and CLK k can 
be estimated by calculating a correlation function between the timing jitter 
5 sequences A^fn] and A(j) k [n], and by searching an offset value at which the 
correlation function shows the maximum value. In summary, the amount of 
offset n 0 ff se t of the clock edges can be obtained from an offset position at which 
the cross-correlation between the instantaneous phase noises show the maximum 
value. 

10 In addition, a deterministic clock skew value can also be obtained by 

obtaining zero-crossing time points of respective clock signals under 
measurement, and by calculating a mean value of time differences between 
corresponding zero-crossings. 

Therefore, according to one mode of the clock skew measurement 

15 method of the present invention, timing jitters A^fn] and A(j) k [n] of the respective 
two clock signals under measurement xj(t) and x k (t) shown in Fig. 6 are first 
obtained, and then a deterministic clock skew value Ts ke J' k between the two 
clock signals under measurement xj(t) and x k (t) is obtained. Next, a timing 
difference between clock edges is obtained by calculating a difference between 

20 the timing jitter sequences A^fn] and A<|) k [n], and then the random component of 
the clock skew between the clock signals under measurement xj(t) and x k (t) is 
obtained. A clock skew Tskew^tn] between the clock signals under 
measurement is obtained by obtaining a sum of the random component of the 
clock skew and the deterministic component T Skew j ' k An obtained clock skew 

25 T skew j,k [n] is shown in Fig. 7. If necessary, an RMS value and a peak-to-peak 
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value of clock skew are obtained from the clock skew sequence T S kew ' k [n] . An 
RMS value TskewjW* of clock skew is a standard deviation of the clock skew 
Tskew^M, and is obtained by the following equation. 



T Skew,RMS 



'n=l 



[sec] 



(22) 



10 



In this case, N is the number of samples of the measured clock skew data, and 
T'skew j ' k is a mean value. In addition, a peak-to-peak value TskewV* °f clock 
skew is a difference between the maximum value and the minimum value of the 
Tskew *[n], and can be obtained by the following equation. 



TiL, P P=max 0 (T4LM)-min n (Tj i L[k]) [sec] (23) 

15 Fig. 8 shows a histogram of clock skew measured by this clock skew 
measurement method. 

According to another aspect of the clock skew measurement method of 
the present invention, a clock skew between clock signals each having a different 
frequency from one another can also be measured. Now, in a clock distribution 

20 network shown in Fig. 9, a system clock signal CLKg from an external system 
clock source 14 is inputted to a clock signal source 11 comprised of a PLL 
(Phase Locked Loop) circuit, where its frequency is multiplied by M. A clock 
signal CLKg whose frequency is M multiple of that of the CLKg is distributed as 
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the clock signals CLKj and CLKk to a network, for example, the registers 12j and 
12k. Fig. 10a shows the system clock signal CLKg, Fig. 10b shows an ideal 
clock signal that is created by multiplying the frequency of the CLKg by M, and 
Fig. 10c shows the distributed clock signal CLKj that is created by multiplying 
5 the frequency of the CLKg by M. A®[1] [rad] of the system clock signal CLK G 
represents a timing fluctuation of its edge from an ideal clock edge. Therefore, 
the clock shown in Fig. 10b that is created by multiplying the frequency of the 
system clock signal CLK G shown in Fig. 10a by M is, as shown in Fig. 10b, a 
clock of a case in which the system clock signal CLKg is multiplied by 2 in the 

10 case of illustrated example. Therefore, the number of clock edges is doubled in 
the case of Fig. 10b. In this case, it is sufficient for a rising edge of a newly 
increased clock to assign the jitter A0[1] of a rising edge of the original system 
clock. In the case where the frequency is multiplied by M, one to one 
relationship of A0[{n/M}] and A^n] can be obtained when the A0[1] is 

15 assigned (M-l) times. In this case, {x} represents the maximum integer that 
does not exceed x. When a clock skew between the clock signals CLKj and 
CLKg is obtained using the equation (6), an equation (24) is obtained. 



T G 'j rnl = T G 'j + 
1 Skewl n J T Skew + 



A<|> j [n] 



-A0 



n 



2tu 



[sec] (24) 



20 



A deterministic clock skew value iskew between the clock signals CLKj and 
CLKg is expressed by a time difference between an ideal clock edge (nMT) j of 
the clock signal CLKj and an ideal clock edge (nMT) g of the system clock 
signal CLK G , and can be obtained by the following equation from the initial 
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phase angles of the respective clock signals. 



-<£H£M^MS tseci (25) 



5 In this case, since the clock signal CLKj is a clock signal that is created by 

multiplying the frequency of the system clock signal CLKg by M, a fundamental 
period T G of the clock signal CLKg is equal to M times of a fundamental period 
Tj of the clock CLKj (T G =MTj). 

In addition, according to further another aspect of the clock skew 

1 0 measurement method of the present invention, a clock skew between the clock 
signals CLKj and CLKk can be measured, using an apparatus that can 
simultaneously measure two channels, by simultaneously sampling only the 
clock signals CLKj and CLK g at first, and then by simultaneously sampling only 
the clock signals CLK k and CLK g . 

15 That is, only the clock signals CLKj and CLK g are simultaneously 

sampled at first, and then a clock skew between the clock signals CLK, and 
CLKg is obtained by the following equation (26) using the second term in the 
right side of the first equation in the equation (6). 
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Next, only the clock signals CLKk and CLK g are simultaneously sampled 
at first, and then a clock skew between the clock signals CLKk and CLK g is 
5 similarly obtained by the following equation (27) 



Tskew M ~ T Skew + 



[n] -S. - A ^ [n] 

\1%J 



(T } 



K 2%j 



[sec] (27) 
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Finally, a difference between the clock skew sequences T S kew 8,i and Tskew 8 ' k 
obtained by the equations (26) and (27), respectively is calculated to estimate a 
clock skew between the clock signals CLKj and CLKk using the following 
equation (28). 



15 



-23- 



T 4ewW = T S^LM- T S^wW 
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(28) 



As a result, when clock skews among N clock signals are measured, the 
required number of simultaneous samplings is decreased to (N-l) two- channel 
5 simultaneous measurements from nC 2 (=N(N-l)/2). In addition, this method 
requires only the minimum number of pins when distributed clock signals on a 
semiconductor chip, for example, are taken out to the outside of the chip. 
Therefore, this method is optimum for an evaluation test of a VLSI. In addition, 
the above procedure can be also applied to clock signals each having different 
10 frequency from one another. 

The clock skew measurement method according to the present invention 
can be used not only for the estimation of a clock skew between distributed 
clocks in a microprocessor unit as described above, but can also be applied to the 
estimation of a clock skew of other signals. 



15 
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TDVUNG JITTER ESTIMATION METHOD 

Next, a timing jitter estimation method used in the clock skew 
measurement method according to the present invention will be explained. 

A jitter-free clock signal is a square wave having a fundamental 
5 frequency fo. This signal can be decomposed by Fourier analysis to harmonics 
having frequencies f 0 , 3f 0 , 5f 0 , — . Since a jitter corresponds to a fluctuation of a 
fundamental frequency of a clock signal under measurement, only signal 
components around the fundamental frequency are handled in the jitter analysis. 

A fundamental component of a jittery clock signal (clock signal under 
10 measurement) can be expressed as a sinusoidal wave, assuming that its 
amplitude is A and its fundamental period is T 0 , by the following equation. 



x(t) = Acos((|>(t)) = A cos 



t + (|> 0 -A<|>(t) 



(29) 



1 5 Here, (|)(t) is an instantaneous phase of the clock signal under measurement, and 
can be expressed by a sum of a linear instantaneous phase component 27rt/T 0 
containing the fundamental period T 0 , an initial phase angle (can be made zero 
in the calculation), and an instantaneous phase noise component A<j)(t). 

When the instantaneous phase noise component A<|)(t) is zero, an interval 

20 between rising zero-crossing points of the clock signal under measurement is a 
constant period T 0 . The non-zero A(J>(t) fluctuates the zero-crossing points of 
the clock signal under measurement. That is, a A<|) (nTo) at a zero-crossing 
point nTo represents a time fluctuation of a zero-crossing point, and is called a 
jitter. Therefore, an instantaneous phase noise A<|>(t), at zero-crossing points, a 
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timing jitter of the clock signal under measurement can be obtained, by 
estimating an instantaneous phase (|>(t) of the clock signal under measurement, 
and by obtaining a difference between an instantaneous phase and a linear phase 
(corresponds to a phase waveform of a jitter-free ideal clock signal) 27rt/To+<t>o. 

5 In the timing jitter estimation method used in the present invention, for 

example, a clock signal under measurement x(t) is transformed into a complex 
analytic signal z(t) at first, and an instantaneous phase §(t) of the clock signal 
under measurement x(t) is estimated from the analytic signal z(t). A linear line 
fitting by least squares method is applied to the estimated instantaneous phase 

10 waveform data to obtain a linear instantaneous phase <|>iinear(t) corresponding to 
an instantaneous phase waveform of a jitter-free ideal signal. Then a difference 
between the instantaneous phase and the linear instantaneous phase (|)iinear(t) is 
calculated to obtain an instantaneous phase noise A<|>(t) of the clock signal under 
measurement. The instantaneous phase noise waveform A(|)(t) is sampled at 

1 5 timings (approximated zero-crossing points) closest to each zero-crossing point 
of a real part x(t) of the analytic signal z(t) to estimate an instantaneous phase 
noise at a zero-crossing timing nTo, i.e., a timing jitter A<|>[n](=A<|)(nTo)). The 
method of obtaining, in this manner, an instantaneous phase noise A4>(t) to 
estimate a timing jitter A(|)[n] was proposed by the inventors of the present 

20 invention, and is described in, for example, "Extraction of Peak-to-Peak and 
RMS Sinusoidal Jitter Using an Analytic Signal Method" by T.J. Yamaguchi, M. 
Soma, M. Ishida, T. Watanabe, and T. Ohmi, Proceedings of 18th IEEE VLSI 
Test Symposium, pp. 395-402, 2000. 

In this timing jitter estimation method, a timing jitter can be estimated 

25 with high accuracy by removing, using waveform clipping means, amplitude 
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modulation (AM) components of the clock signal under measurement in the state 
that phase modulation components corresponding to jitter of the clock signal 
under measurement are kept in the clock signal under measurement. In 
addition, it is desirable to remove low frequency components of the 
5 instantaneous phase noise using low frequency component removing means. 

INSTANTANEOUS PHASE ESTIMATION METHOD USING ANALYTIC 
SIGNAL 

An analytic signal z(t) of a clock signal under measurement x(t) is 
10 defined by a complex signal expressed by the following equation (30). 

z(t)-x(t)+ji(t) (30) 

In this case, j represents an imaginary unit, and an imaginary part x (t) of the 
1 5 complex signal z(t) is a Hilbert transform of a real part x(t). 

On the other hand, Hilbert transform of a time waveform x(t) is defined 
by the following equation (31). 

x(t) = H[x(t)] = Ij^dt (31 ) 

71 J ~ G0 t-T 

20 

In this case, x (t) is a convolution of the function x(t) and (l/7if). That is, 
Hilbert transform is equivalent to an output at the time when the x(t) is passed 
through a all pass filter. However, the output x (t) in this case has not been 
changed in terms of its spectrum components, but its phase has been shifted by 
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7l/2. 

Analytic signal and Hilbert transform are described in, for example, 
"Probability, Random Variables, and Stochastic Processes" by A. Papoulis, 2nd 
edition, McGraw-Hill Book Company, 1984. 
5 An instantaneous phase waveform (|)(t) of a clock signal under 

measurement x(t) can be obtained from an analytic signal z(t) using the 
following equation (32). 



10 



(|>(t) = tan 



-1 



x(t) 
x(t) 



(32) 



Next, an algorithm for estimating an instantaneous phase using Hilbert 
transform will be described. First, Hilbert transform is applied to a clock signal 
under measurement x(t). 
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x(t) = A cos 



rr 2% \ 



t+* 0 -A*(t) 



(33) 



Then a signal x(t) corresponding to an imaginary part of an analytic signal z(t) 
is obtained as follows. 



20 



x(t) = H[x(t)] = Asin 



"2^ 
W T oj 



t+4> 0 -A4>(0 



(34) 



Then the clock signal under measurement x(t) is transformed into an analytic 
signal z(t) as follows. 
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z(t) = x(t) + jx(t) = A cos 



W T 0J 



t + <|) 0 -A(j)(t) 



+ jAsin 



W T oj 



t+<|>o-A<Kt) 



(35) 
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Here, a band-pass filtering process has been applied to the obtained analytic 
signal. This is because a jitter corresponds to a fluctuation of a fundamental 
frequency of a clock signal under measurement and hence only signal 
components around the fundamental frequency are processed in a jitter analysis. 
Next, a phase function <|)(t) is estimated from the obtained analytic signal z(t) 
using the equation (32). 



<Kt) = 



^t+4> 0 -A4>(t) 



mod 2% 



[rad] 



(36) 



Here, the <|)(t) is expressed using principal values of phase in the range of -% to 
+7C, and has a discontinuity point at the proximity of a point where the phase 
1 5 changes from +n to -it. Finally, by unwrapping (that is, integer multiples of 2n 
are appropriately added to principal values (])(t)) the discontinuous phase function 
<Kt), a continuous instantaneous phase function §(t), from which discontinuity 
has been removed, can be obtained. 



20 



<Kt) = 



V 1 0J 



t+<j> 0 -A<|>(t) 



[rad] 



(37) 
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A phase unwrapping method is described in "The Cephstrum: A Guide to 
Processing" by Donald G. Childers, David P. Skinner, and Robert C. Kemerait, 
Proceedings of IEEE, vol. 65, pp. 1482-1442, 1977. 

5 TRANSFORMATION TO ANALYTIC SIGNAL USING FAST FOURIER 
TRANSFORM 

A transformation from a clock signal under measurement to an analytic 
signal can be realized by a digital signal processing using a transformation of a 
signal in time domain into a signal in frequency domain like Fast Fourier 

10 Transformation. 

First, FFT is applied to a digitized clock signal under measurement x(t) 
shown in Fig. 11 to obtain a both-sided spectrum (having positive and negative 
frequencies) X(f). Fig. 12A shows the obtained both-sided spectrum X(f). 
Next, as shown in Fig. 12B, only data around the fundamental frequency in the 

1 5 positive frequency components of the spectrum X(f) are retained, and the 
remaining data are replaced by zeros. In addition, the positive frequency 
components are doubled. These processes in frequency domain correspond to 
limiting bandwidth of the clock signal under measurement and transforming the 
clock signal under measurement into an analytic signal Z(f) in time domain. 

20 Finally, by applying inverse FFT to the band-limited signal Z(f), a band-limited 
analytic signal z(t) can be obtained. 

Transformation to an analytic signal using FFT is described in, for 
example, "Random Data: Analysis and Measurement Procedure" by 
J.S. Bendat and A.G. Piersol, 2nd edition, John Wiley & Sons, Inc., 1986. 

25 In addition, when an estimation of instantaneous phase is the object of 
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the process, the process of doubling the positive frequency components can be 
omitted. 

APPROXIMATED ZERO-CROSSING POINT DETECTION METHOD 
Next, an approximated zero-crossing point detection method will be 

5 described. First of all, the maximum value of a real part x(t) of an analytic 
signal of the inputted clock signal under measurement is defined as 100% level, 
and the minimum value is defined as 0% level to calculate 50% level signal 
value V5o% as a zero crossing level. A difference between a sample value and 
the 50% level Vso% and a difference between its adjacent sample value and the 

10 50 % level V 50 % of the x(t), i.e., (x(j-l)-V 50 %) and ( x(j)-V 50 o /o ) are calculated, and 
furthermore a product of those difference values (x(j-l)-V5o%) x ( x(j)-V 50 %) is 
calculated. When the x(t) crosses 50% level, i.e., zero-crossing level, the sign 
of its sample value (x(j-l) -V 50 %) or (xG)-V 50 %) changes from a negative sign to a 
positive sign or from a positive sign to a negative sign. Therefore, when the 

1 5 product is negative, it is detected that the x(t) has passed the zero-crossing level, 
and a time point j-1 or j at which a smaller absolute value of the sample value 
(x(j-l) -V5o%) or (x(j)-V5oo/ 0 ) is detected is obtained as an approximated zero 
crossing point. 

20 WAVEFORM CLIPPING 

Waveform clipping means removes AM components from an input signal 
and retains only PM components in the input signal. A waveform clipping is 
performed by: 1) multiplying an analog or digital input signal by a constant, 2) 
replacing a signal value greater than a predetermined threshold value Thl with 

25 the threshold value Thl, and 3) replacing a signal value less than a 
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predeteraiined threshold value Th2 with the threshold value Th2. Here, it is 
assumed that the threshold value Thl is greater than the threshold value Th2. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a diagram showing an example of a clock skew measurement 
using a time interval analyzer; 

Fig. 2 is a diagram typically showing a clock distribution network; 

Fig. 3 is a diagram typically showing a timing of clock skew; 

Fig. 4 is a diagram typically showing a relationship between a timing 
jitter and a clock skew; 

Fig. 5 A is a diagram showing an example of a timing jitter A$[n] of a 
clock signal under measurement xj(t); 

Fig. 5B is a diagram showing an example of a timing jitter A<|) k [n] of a 
clock signal under measurement Xk(t); 

Fig. 6 is a diagram showing an example of a clock signal under 
measurement; 

Fig. 7 is a diagram showing an example of a clock skew between clock 
signals under measurement measured by a clock skew measurement method 
according to the present invention; 

Fig. 8 is a diagram showing an example of a histogram of clock skew 
between clock signals under measurement measured by the clock skew 
measurement method according to the present invention; 

Fig. 9 is a diagram typically showing a clock distribution network having 
different clock domains; 

Fig. 10 is a diagram typically showing a principle of clock skew 



-32- 



measurement using a frequency multiplication; 

Fig. 11 is a diagram showing an example of a digitized clock signal under 
measurement; 

Fig. 12 A is a diagram showing an example of a both-sided power 
5 spectrum of a clock signal under measurement obtained by FFT; 

Fig. 12B is a diagram showing an example of a band-limited single-sided 
power spectrum; 

Fig. 13 is a block diagram showing a functional configuration of an 
embodiment of a clock skew measurement apparatus according to the present 
10 invention; 

Fig. 14 is a flow-chart showing an embodiment of a clock skew 
measurement method according to the present invention; 

Fig. 15 is a block diagram showing a functional configuration of a 
specific example of a deterministic clock skew estimator 102 shown in Fig. 13; 
15 Fig. 16 is a flow-chart showing a processing example of a deterministic 

clock skew estimating step of the step 202 shown in Fig. 14; 

Fig. 17 is a block diagram showing a functional configuration of another 
embodiment of the clock skew measurement apparatus according to the present 
invention; 

20 Fig. 18 is a flow-chart showing another embodiment of the clock skew 

measurement method according to the present invention; 

Fig. 19 is a block diagram showing an example of a functional 

configuration of a timing jitter estimator used in the clock skew measurement 

apparatus according to the present invention; 
25 Fig. 20 is a flow-chart showing an example of a timing jitter estimation 
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method used in the clock skew measurement method according to the present 
invention; 

Fig. 21 is a block diagram showing another example of a functional 
configuration of an analytic signal transformer used in the clock skew 
5 measurement apparatus according to the present invention; 

Fig. 22 is a block diagram showing further another example of a 
functional configuration of an analytic signal transformer used in the clock skew 
measurement apparatus according to the present invention; 

Fig. 23 is a flow-chart showing further another example of the analytic 
10 signal transformation method used in the clock skew measurement method 
according to the present invention; 

Fig. 24 is a block diagram showing a portion of another embodiment of 
the apparatus according to the present invention; 

Fig. 25 is a flow-chart showing a portion of another embodiment of the 
1 5 method according to the present invention; and 

Fig. 26 is a block diagram showing another example of a functional 
configuration of the deterministic clock skew estimator 102; 

DESCRIPTION OF THE EMBODIMENTS 

20 Embodiments of the present invention are described below. 

Fig. 13 shows a functional configuration of an embodiment of a clock 
skew measurement apparatus according to the present invention. This clock 
skew measurement apparatus 100 comprises timing jitter estimators 101a and 
101b for estimating timing jitter sequences A^fn] and A<j) k [n] of clock signals 

25 under measurement xj(t) and x k (t), respectively, a deterministic clock skew 



-34- 

estimator 102 for estimating a timing error between ideal clock edges of the 
respective clock signals under measurement and for estimating a deterministic 
component Taw* of clock skew, a clock skew estimator 103 to which the timing 
jitter sequences A^n] and A<|) k [n] are inputted for calculating a timing difference 

5 sequence between those timing jitter sequences A(^[n] and A(|> k [n] to output a 
clock skew sequence Tskew^M, and a clock skew detector 104 for obtaining a 
clock skew between the clock signals under measurement from the clock skew 
sequence. In addition, the clock skew detector 104 comprises a peak-to-peak 
detector 105 for obtaining a difference between the maximum value and the 

10 minimum value of the clock skew sequence T S kew [n], an RMS detector 106 for 
calculating an RMS value as a standard deviation of the clock skew sequence, 
and a histogram estimator 107 for obtaining a histogram of the clock skew 
sequence. The timing jitter estimators 101a and 101b estimate initial phase 
angles §j and <|>o k of the clock signals under measurement Xj(t) and Xk(t), 

1 5 respectively in addition to the timing j itter sequences A^ [n] and A<|) k [n] , and 
output the estimated initial phase angles to the deterministic clock skew 
estimator 102. A specific configuration of the timing jitter estimators 101a and 
101b will be described later. 

Next, the operation in the case where a clock skew between the clock 

20 signals under measurement xj(t) and x k (t) is measured using the clock skew 
measurement apparatus 100 of this embodiment will be described. Fig. 14 
shows a processing procedure of the embodiment of the clock skew 
measurement method according to the present invention. First, in step 201 , 
initial phase angles fyj and (j) 0 k , and timing jitter sequences A^fn] and A<|) k [n] of 

25 the respective clock signals under measurement Xj(t) and Xk(t) are estimated by 
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the timing jitter estimators 101a and 101b, respectively. Next, in step 202, a 
difference between the initial phase angles V and <|>o k of the respective clock 
signals under measurement is calculated by the deterministic clock skew 
estimator 102 to estimate a deterministic component TskeJ* of clock skew 
5 between the clock signals under measurement. Next, in step 203, a clock skew 
sequence T S kew j,k |n] between the clock signals under measurement xj(t) and x^t) 
is estimated by the clock skew estimator 103 from the timing jitter sequences 
A^Jn] and A<|) k [n], and the deterministic component xskew' k of clock skew. 
Finally, in step 204, a clock skew value between the clock signals under 
1 0 measurement xj(t) and x k (t) is obtained by the clock skew detector 1 04 from the 
estimated clock skew sequence T S kew' k , and the process ends. 

In the step 202 for estimating a deterministic component of clock skew 
between the clock signals under measurement, the deterministic clock skew 
estimator 102 obtains a deterministic component of clock skew between the 
1 5 clock signals under measurement using the equation (16). In addition, in the 
step 202, the deterministic clock skew estimator 102 may obtain, if necessary, an 
absolute value of the equation (16). In addition, in the step 203 for estimating a 
clock skew sequence between the clock signals under measurement, the clock 
skew estimator 103 obtains a clock skew sequence T sk ew j ' k [n] between the clock 
20 signals under measurement using the equation (6). In the step 204 for obtaining 
a clock skew value between the clock signals under measurement, the 
peak-to-peak detector 105 obtains a peak-to-peak value of clock skew using the 
equation (23), the RMS detector 106 obtains an RMS value of clock skew using 
the equation (22), and the histogram estimator 107 obtains a histogram from the 
25 clock skew sequence. Alternatively, an RMS value and/or the peak-to-peak 
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value may be obtained from only the second term of the equation (6). In 
addition, the step 201 for estimating initial phase angles and timing jitter 
sequences of the clock signals under measurement may be replaced by the 
processing procedure shown in Fig. 20. Furthermore, the step 202 for 
estimating a deterministic component of clock skew between the clock signals 
under measurement may be replaced by the processing procedure shown in Fig. 
16. 

The clock skew measurement apparatus shown in Fig. 13 can also be 
modified as an apparatus for estimating only random component of clock skew. 
In this case, the deterministic clock skew estimator 102 for obtaining a 
deterministic component of clock skew is omitted. Similarly, the clock skew 
measurement method shown in Fig. 14 can also be modified as a method of 
estimating only random component of clock skew. In this case, the step 202 for 
estimating a deterministic component of clock skew from initial phase angles of 
the clock signals under measurement is omitted. 

The deterministic clock skew estimator 102 shown in Fig. 13 estimates a 
deterministic component of clock skew from a difference between initial phase 
angles <j) 0 j and <|)o k of the respective clock signals under measurement. However, 
the deterministic clock skew estimator can also be materialized by the 
configuration shown in Fig. 15. That is, this deterministic clock skew estimator 
102 to which initial phase angles <J> 0 j and (|>o k and timing jitter sequences A#[n] 
and A(j) k [n] of the respective clock signals under measurement xj(t) and x k (t) are 
inputted comprises an offset estimator 301 for estimating an offset n 0 ff S et between 
corresponding clock edges of the respective clock signals under measurement 
from those timing jitter sequences A^tn] and A<l> k [n], and a deterministic clock 
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skew calculator 302 for calculating a deterministic component Tsew* of clock 
skew between the clock signals under measurement from the initial phase angles 
V and <|>o k and the offset noff se t between the corresponding clock edges estimated 
by the offset estimator 301. 

5 The operation in the case where a deterministic component of clock skew 

between the clock signals under measurement is estimated using this 
deterministic clock skew estimator 102 will be described. Fig. 16 shows the 
processing procedure. First, in step 40 1 , from the inputted timing jitter 
sequences A$[n] and A<|> k [n] of the clock signals under measurement, an offset 

1 0 position, at which correlation function between those timing jitter sequences 
shows the largest value, is obtained by the offset estimator 301 to estimate an 
offset noffset between the corresponding clock edges. Next, in step 402, a 
deterministic component T S kew j ' k of clock skew between the clock signals under 
measurement xj(t) and x k (t) is calculated by the deterministic clock skew 

1 5 calculator 302 from the inputted initial phase angles V and <|>o k and the offset 
n 0 ffset between the corresponding clock edges, and then the process ends. In the 
step 402 for calculating a deterministic component of clock skew between the 
clock signals under measurement, the deterministic clock skew calculator 302 
obtains a deterministic component of clock skew between the clock signals 

20 under measurement using the equation (21). As indicated by dashed lines in 
Fig. 15, instantaneous phase noises A^t) and A<|) k (t) to be explained in Fig. 19 
may be inputted to the offset estimator 301 to obtain an offset position at which a 
correlative value between those A^t) and A<|> k (t) becomes the largest value for 
estimating an offset no ffse t. 

25 Fig. 17 shows a functional configuration of another embodiment of the 
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clock skew measurement apparatus according to the present invention. This 
clock skew measurement apparatus 500 comprises timing jitter estimators 101a, 
101b, 101c, and lOld for respectively estimating timing jitter sequences A^fn], 
A(|) 8 [n], A()) k [n], and A<|) s [n] of respective clock signals under measurement Xj(t), 

5 x g (t), Xk(t), and x g (t), deterministic clock skew estimators 102a and 102b for 
respectively estimating a timing error Et g,i between ideal clock edges of the clock 
signals under measurement Xj(t) and x g (t) and a timing error Et gJ between ideal 
clock edges of the clock signals under measurement Xk(t) and x g (t) to 
respectively estimate deterministic components Tskew 8 "* and T S kew 8 ' k using those 

10 timing errors Et gJ and Et g ' k , clock skew estimators 103a and 103b to which the 
timing jitter sequences A^n], A<|> g [n], and A<|) k [n], A(|> g [n] are respectively 
inputted for calculating timing difference sequences between the two inputs at 
each clock skew estimator respectively to output clock skew sequences T skew 8 [n] 
and Tskew^W respectively, a clock skew estimator 501 to which those clock 

1 5 skew sequences T skew 8j [n] and T ske w 8 ' k [n] are inputted for obtaining a difference 
between those clock skew sequences to estimate a clock skew sequence 
Tskew^W, and a clock skew detector 104 for obtaining a clock skew value 
between the clock signals under measurement from the clock skew sequence 
obtained by the clock skew estimator 501 . For simplicity, the explanation of 

20 portions duplicated with those in Fig. 13 is omitted. 

Next, the operation in the case where a clock skew between the clock 
signals under measurement is measured using the clock skew measurement 
apparatus 500 according to the present invention will be described. Fig. 18 
shows a processing procedure of the clock skew measurement method according 

25 to the present invention. First, in step 601 , initial phase angles (|>o J and (j) 0 g , and 
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timing jitter sequences A^n] and A(|) g [n] of the respective clock signals under 
measurement xj(t) and Xg(t) are estimated by the timing jitter estimators 101a and 
101b, respectively. Next, in step 602, a difference between the initial phase 
angles <j) 0 j and (J) 0 8 of the clock signals under measurement is calculated by the 
5 deterministic clock skew estimator 1 02a to estimate a deterministic component 
fskew 8j of clock skew between the clock signals under measurement. Next, in 
step 603, a clock skew sequence T skew gJ [n] between the clock signals under 
measurement is estimated by the clock skew estimator 103a from the timing 
jitter sequences A^n] and A<|) 8 [n], and the deterministic component x S kew 8J of 
10 clock skew. 

Next, in step 604, initial phase angles <j> 0 k and (j) 0 8 , and timing jitter 
sequences A<|) k [n] and A<|) 8 [n] of the respective clock signals under measurement 
x k (t) and x g (t) are estimated by the timing jitter estimators 101c and lOld, 
respectively. Next, in step 605, a difference between the initial phase angles ty* 
15 and (J) 0 8 of the clock signal under measurement obtained in the step 604 is 
calculated by the deterministic clock skew estimator 102b to estimate a 
deterministic component x S kew s * of clock skew between the clock signals under 
measurement. Next, in step 606, a clock skew sequence T S kew g,k [n] between the 
clock signals under measurement is estimated by the clock skew estimator 103b 
20 from the timing jitter sequences A(|> k [n] and A(j) 8 [n] obtained in the step 604 and 
the deterministic component T S kew 8 ' k of clock skew obtained in the step 605. 
Next, in step 607, a clock skew sequence T S kew' k [n] between the clock signals 
under measurement Xj(t) and x k (t) is estimated by the clock skew estimator 501 
from the clock skew sequences T S kew 8J [n] and T skew 8 ' k [n] respectively obtained in 
25 the steps 603 and 606 . Finally, in step 608, a clock skew value between the 
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clock signals under measurement xj(t) and x k (t) is obtained by the clock skew 
detector 104 from the clock skew sequence T S kew J,k [n] estimated in the step 607, 
and the process ends. In the step 607 for estimating the clock skew sequence 
between the clock signals under measurement xj(t) and x k (t), the clock skew 
5 estimator 501 obtains a clock skew sequence between the clock signals under 
measurement using the equation (28). The process sequence of the steps 
60 1-603 and the steps 604-606 may be exchanged. For simplicity, the 
explanation of portions duplicated with those in Fig. 14 is omitted. 

The clock skew measurement apparatus shown in Fig. 17 may also be 
1 0 configured as an apparatus for estimating only random component of clock skew. 
In this case, the deterministic clock skew estimators 102a and 102b for obtaining 
the deterministic components of clock skew are omitted. Similarly, the clock 
skew measurement method shown in Fig. 18 may also be modified as a method 
of estimating only random component of clock skew. In this case, the steps 602 
15 and 605 for estimating deterministic components of clock skew from the initial 
phase angles of the clock signals under measurement are omitted. 

As indicated by dashed lines in Fig. 13, a timing jitter sequence estimated 
by one of the timing jitter estimators 101a and 101b, i.e., the timing jitter 
sequence A<|) k [n] estimated by the timing jitter estimator 101b in the case of Fig. 
20 13 may be copied (M-l) times by the frequency multiplier 701 to obtain a timing 
jitter sequence that is to be obtained when the frequency of the clock signal 
under measurement x k (t) is multiplied by M, and to supply the obtained timing 
jitter sequence to the clock skew estimator 103. By such a process, in the clock 
distribution system previously described with reference to Fig. 9, the clock 
25 signal x k (t) in Fig. 13 corresponds to the system clock signal CLKg in Fig.9, and 
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the clock signal xj(t) is a signal created by multiplying the frequency of the clock 
signal CLK G (x k (t)) by M. In this case, a clock skew T S kew Gj [n] between the 
clock signal CLKo^t)) and the clock signal xj(t) can be obtained by the 
equation (24). 

5 Regarding the processing procedure of the clock skew measurement 

method in this case, as shown in Fig. 14, after timing jitter sequences are 
obtained in the step 201, the timing jitter sequence estimated by the timing jitter 
estimator 101b may be assigned, for example (M-l) times as indicated by dashed 
lines, in step 801 by the frequency multiplier to obtain a timing jitter sequence 

10 that is to be obtained when the frequency of the clock signal under measurement 
is multiplied by M, and then the process may move to the step 202. At this 
time, in the step 202 for estimating a deterministic component of clock skew 
between the clock signals under measurement, the deterministic clock skew 
estimator 102 obtains a deterministic component of clock skew between the 

1 5 clock signals under measurement using the equation (25). In addition, in the 
step 203 for estimating a clock skew sequence between the clock signals under 
measurement, the clock skew estimator 103 obtains a clock skew sequence 
between the clock signals under measurement using the equation (24). 

Also in the case of using the frequency multiplier 701, when the clock 

20 skew measurement apparatus is constructed as an apparatus for estimating only 
random component of clock skew, the deterministic clock skew estimator 102 for 
obtaining a deterministic component of clock skew may be omitted. In this 
case, in the clock skew measurement method, the step 202 for estimating a 
deterministic component of clock skew from the initial phase angles of the clock 

25 signals under measurement may be omitted. 
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In addition, the frequency multiplier 701 may also be built in the clock 
skew measurement apparatus shown in Fig. 17. In this case, the frequency 
multiplier is inserted in series in each output side of the timing jitter estimators 
101b and 10 Id. Similarly, a step of multiplying a frequency can also be added 

5 to the clock skew measurement method shown in Fig. 18. In this case, the step 
of multiplying a frequency is inserted after each of the steps 601 and 604 for 
estimating timing jitters. 

Fig. 19 shows an example of configuration of each of the timing jitter 
estimators 101a, 101b, 101c, and lOld. This timing jitter estimator 900 is 

10 described in, for example, "Extraction of Peak-to-Peak and RMS Sinusoidal 
Jitter Using an Analytic Signal Method" by T.J. Yamaguchi, M. Soma, M. 
Ishida, T. Watanabe, and T. Ohmi, Proceedings of 18th IEEE VLSI Test 
Symposium, pp. 395-402, 2000. This timing jitter estimator 900 comprises an 
analytic signal transformer 901 for transforming a clock signal under 

1 5 measurement into a band-limited complex analytic signal, an instantaneous 
phase estimator 902 for obtaining an instantaneous phase of the analytic signal 
transformed by the analytic signal transformer 901, a linear phase remover 903 
for removing a linear instantaneous phase from the instantaneous phase 
estimated by the instantaneous phase estimator 902 to obtain an instantaneous 

20 phase noise, a zero-crossing detector 904 to which a real part of the analytic 
signal is inputted from the analytic signal transformer 901 for generating 
sampling pulses at timings (approximated zero-crossing points) close to 
zero-crossing timings of the real part of the analytic signal, a zero-crossing 
sampler 905 to which the instantaneous phase noise estimated by the linear phase 

25 remover 903 is inputted for sampling the instantaneous phase noise using the 
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sampling pulses from the zero-crossing detector 905. The analytic signal 
transformer 901 may be constructed such that a pass bandwidth of a signal can 
arbitrarily be changed. In addition, the linear phase remover 903 obtains an 
initial phase angle of the clock signal under measurement simultaneously with 

5 the instantaneous phase noise, and outputs the obtained initial phase angle to the 
deterministic clock skew estimator. 

A processing procedure in this timing jitter estimator 900 will be 
explained with reference to Fig. 20. In step 1001, an inputted clock signal 
under measurement is transformed by the analytic signal transformer 901 into an 

10 analytic signal whose predetermined frequency components are selectively 
passed. In step 1002, an instantaneous phase of the clock signal under 
measurement is estimated by the instantaneous phase estimator 902 using the 
analytic signal. In step 1003, a linear instantaneous phase corresponding to an 
ideal clock signal is estimated by the linear phase remover 903 from the 

1 5 instantaneous phase to obtain an initial phase angle of the clock signal under 
measurement, hi step 1 004, the linear instantaneous phase is removed by the 
linear phase remover 903 from the instantaneous phase to estimate an 
instantaneous phase noise A^t). At the same time, in step 1005, timings 
(approximated zero-crossing points) closest to zero-crossing points of a real part 

20 of the analytic signal are detected by the zero-crossing detector 904 from the real 
part of the analytic signal using the previously explained zero-crossing point 
detection method. Finally, in step 1006, only the instantaneous phase noise 
data at the approximated zero-crossing points from the instantaneous phase noise 
are sampled by the zero-crossing sampler 904 to estimate a timing jitter sequence 

25 A^fn], and the process ends. 
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In the analytic signal transformer 901 used in the timing jitter estimator 
900, for example as shown in Fig. 19, only components around a fundamental 
frequency are extracted by a band-pass filter 1101 from the clock signal under 
measurement to band-limit the clock signal under measurement. In addition, 
5 the band-limited clock signal under measurement is inputted to a Hilbert 

transformer 1102 to Hilbert-transform this signal, and an output of the band-pass 
filter 1 1 0 1 is outputted as a real part of the analytic signal. An output of the 
Hilbert transformer 1102 is outputted as a constant part of the analytic signal. 
The band-pass filter 1101 may be either an analog filter or a digital filter, or may 
10 be implemented using a digital signal process such as FFT or the like. In 
addition, the band-pass filter 1101 may be constructed such that the pass 
bandwidth of the signal can arbitrarily be changed. 

Fig. 21 shows another configuration example of the analytic signal 
apparatus 90 1 used in the timing jitter estimator 900. For example, FFT (Fast 
1 5 Fourier Transform) is applied to the clock signal under measurement by a time 
domain to frequency domain transformer 1301 to transform a signal in time 
domain into a both-sided spectrum signal (for example, Fig. 11) in frequency 
domain. Negative frequency components of the transformed both-sided 
spectrum signal in frequency domain are replaced by zero by the bandwidth 
20 limiter 1302 to obtain a single-sided spectrum signal. In addition, regarding, 
this single sided, this single-sided spectrum signal, only components around the 
fundamental frequency of the clock signal under measurement are retained and 
the other frequency components are replaced by zeros to band-limit the signal in 
frequency domain. Inverse FFT is applied by a frequency domain to time 
25 domain transformer 1 303 to the band-limited single-sided spectrum signal to 
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transform the signal in frequency domain into an analytic signal in time domain. 

Fig. 22 shows further another configuration example of the analytic 
signal transformer 90 1 used in the timing jitter estimator 900 . This analytic 
signal transformer 1500 comprises a buffer memory 1501 for storing therein a 
5 clock signal under measurement, a signal extractor 1 502 for extracting the signal 
in the sequential order from the buffer memory 1501 such that the signal being 
extracted is partially overlapped with the signal extracted just before, a window 
function multiplier 1503 for multiplying each extracted partial signal by a 
window function, a time domain to frequency domain transformer 1504 for 
10 transforming each partial signal multiplied by the window function into a 

both-sided spectrum signal in frequency domain, a bandwidth limiter 1505 for 
extracting only components around a positive fundamental frequency of the 
clock signal under measurement from the both-sided spectrum signal 
transformed in frequency domain, a frequency domain to time domain 
1 5 transformer 1 506 for inverse-transforming an output of the bandwidth limiter 
1505 into a signal in time domain, and an inverse window function multiplier 
1507 for multiplying the transformed signal in time domain by an inverse 
number of the window function to obtain a band-limited analytic signal. The 
time domain to frequency domain transformer 1504 and the frequency domain to 
20 time domain transformer 1 506 may be packaged using FFT and inverse FFT, 
respectively. In addition, the bandwidth limiter 1505 may be constructed such 
that the pass bandwidth of the signal can arbitrarily be changed. 

The operation in the case of transforming the clock signal under 
measurement into a band-limited analytic signal using this analytic signal 
25 transformer 1 500 will be described with reference to Fig. 23 . First, in step 
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1601, the buffer memory 1501 stores therein the clock signal under 
measurement. Next, in step 1602, the signal extractor 1502 extracts a portion 
of the stored signal from the buffer memory 1501 . In step 1603, the window 
function multiplier 1503 multiplies the extracted sectioned signal by the window 

5 function. In step 1604, the time domain to frequency domain transformer 1 504 
applies FFT to the partial signal multiplied by the window function to transform 
the signal in time domain into a both-sided spectrum signal in frequency domain. 
In step 1605, the bandwidth limiter 1505 replaces negative frequency 
components of the transformed both-sided spectrum signal in frequency domain 

10 with zero. Furthermore, in step 1606, the bandwidth limiter 1505 retains only 
components around the fundamental frequency of the clock signal under 
measurement in the single-sided spectrum signal produced by replacing negative 
frequency components with zero and replaces the other frequency components 
with zero to limit the bandwidth of the signal in frequency domain. In step 

15 1 607, the frequency domain to time domain transformer 1 506 applies inverse 
FFT to the band-limited single-sided spectrum signal in frequency domain to 
transform the signal in frequency domain into a signal in time domain. In step 
1608, the inverse window function multiplier 1507 multiplies the 
inverse-transformed signal in time domain by an inverse number of the window 

20 function used in the multiplication in the step 1603 to obtain a band-limited 
analytic signal. Finally, in step 1609, a check is made to see if there is 
unprocessed data remaining in the buffer memory 1603. If there is unprocessed 
data remaining in the buffer memory, the signal extractor 1502 extracts, in step 
1610, the signal in the sequential order from the buffer memory 1501 such that 

25 the signal being extracted is partially overlapped with the signal extracted just 
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before, and thereafter the steps 1603, 1604, 1605, 1606, 1607, 1608 and 1609 are 
repeated. If there is no unprocessed data, the process ends. The processing 
sequence of the step 1605 and the step 1606 may be exchanged. That is, only 
components around the fundamental frequency of the signal under measurement 
5 are retained and the other frequency components are replaced with zeros first to 
limit the bandwidth of the signal in frequency domain, and thereafter negative 
frequency components of the both-sided spectrum signal may be replaced with 
zero. 

In the linear phase remover 903 in the timing jitter estimator 900 shown 

10 in Fig. 19, for example as illustrated in the figure, an inputted instantaneous 
phase is converted into a continuous instantaneous phase by a continuous phase 
converter 91 . Regarding the continuous instantaneous phase, its instantaneous 
linear phase, i.e., a linear instantaneous phase corresponding to a jitter-free ideal 
signal is estimated by a linear phase estimator 92 using, for example, a linear 

1 5 trend estimation method, i.e., by applying a linear line fitting by least squares 
method to the continuous instantaneous phase, to output an initial phase angle 
(t> 0 j (t) of the clock signal under measurement xj(t). In addition, the linear 
instantaneous phase is subtracted from the continuous instantaneous phase by a 
subtractor 93 to output an instantaneous phase noise A<])(t). 

20 Further, Fig. 19, Fig. 21, and Fig. 22 are shown in the international 

publication WO00/46606 (published on August 10, 2000) report. 

As indicated by dashed lines in Fig. 13, analog clock signals under 
measurement xj(t) and x k (t) may be digitized by AD converters 1701a and 
1701b, respectively to convert the analog clock signals into digital signals, and 

25 the digital signals may be inputted to the respective timing jitter estimators 101a 
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and 10 lb. In addition, as indicated by dashed lines in Fig. 13, waveform 
clippers 1901a and 1901b may be provided to supply those input signals xj(t) and 
x k (t) to the AD converters 1701a and 1701b or the timing jitter estimators 101a 
and 101b, respectively after removing AM components in the state that phase 

5 modulation components, that are jitter components of the respective input 
signals, are retained in the input signals. The waveform clippers 1901a and 
1901b may be provided in the output sides of the AD converters 1701a and 
1701b, respectively. 

In addition, as indicated by dashed lines in Fig. 19, low frequency 

10 components of the instantaneous phase noise may be removed by a low 

frequency component remover 2101 from the instantaneous phase noise A<J)(t) 
outputted from the linear phase remover 903, and the instantaneous phase noise 
from which low frequency components have been removed may be supplied to 
the zero-crossing sampler 905. 

15 In the above description, the instantaneous phase noise A<|)(t) is sampled 

at approximated zero-crossing points to obtain a timing jitter sequence A^n]. 
However, since the linear phase remover 903 has a configuration as shown in 
Fig. 19, for example as indicated by dashed lines in Fig. 24, the sampling at 
approximated zero-crossing points may be inserted in series between the 

20 instantaneous phase estimator 902 and the continuous phase converter 91 . 

Alternatively, the sampling at approximated zero-crossing points may be inserted 
in series between the continuous phase converter 91 and the linear phase 
estimator 92/subtractor 93 . In such a configuration, a timing jitter sequence 
A(^[n] can also be obtained from the subtracter 93. 

25 In addition, since the estimation of an instantaneous phase noise A(j)(t) 
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from an instantaneous phase is performed by the configuration shown by the 
linear phase remover 903 in Fig. 19, the processing procedure is, as shown in 
Fig. 25, that after an instantaneous phase is obtained in the step 1002 in Fig. 20, 
in step 1003a, the instantaneous phase is converted into a continuous 

5 instantaneous phase by the continuous phase converter 91, and in step 1003b, a 
linear instantaneous phase of the continuous instantaneous phase is estimated by 
the linear phase estimator 92 from the continuous instantaneous phase. After 
that, in step 1004, the linear instantaneous phase is removed by the subtracter 93 
from the continuous instantaneous phase to obtain an instantaneous phase noise 

10 A<Kt). 

Therefore, similarly to the configuration shown in Fig. 24, as indicated 
by dashed lines in Fig. 25, the approximated zero-crossing sampling may be 
applied to the instantaneous phase in step 2001 after the step 1002 to obtain a 
sample sequence of the instantaneous phase. And then the process may move 

15 to the step 1003a to transform the sample sequence into a continuous 
instantaneous phase. 

Alternatively, in step 2002, the continuous phase obtained in the step 
1003a may be sampled at the approximated zero-crossing points to obtain a 
sample sequence of the continuous instantaneous phase, and then the process 

20 may move to step 1003b to estimate an instantaneous linear phase from the 

sample sequence of the continuous instantaneous phase. In either case, in step 
1004, there is obtained a timing jitter sequence A$[n] by sampling the 
instantaneous phase noise at the approximated zero-crossing points. 

Regarding the deterministic clock skew estimator 102 in Fig. 14, for 

25 example, as shown in Fig.26, the clock signals under measurement xj(t) and x k (t) 
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may be inputted to zero-crossing timing detectors 81 and 82, respectively. 
Then zero-crossing timing sequences t j zer0 , c (n) and t^er^n) of those respective 
signals may be inputted to a subtracter 83 to obtain a time difference sequence 
between corresponding zero-crossing time points of the zero-crossing timing 

5 sequences tUccO 1 ) 311(1 1?W(n). Th e11 a mean value of those ti™ 6 differences 
may be calculated by a mean value calculator 84 to obtain the mean value as a 
deterministic clock skew value T i ' k S hew. 

The apparatus shown in Fig. 13 and Fig. 17 can be functioned by 
executing programs in a computer. 

10 According to the clock skew measurement apparatus and the clock skew 

measurement method of the present invention, a random spread (random 
components) of clock skew can be measured, and therefore an efficiency of 
clock skew test can greatly be improved. Moreover, if necessary, the apparatus 
can be constructed such that one of the clock skew components that is 

1 5 deterministically determined by clock distribution network paths are also 
simultaneously measured. 

In addition, according to the embodiments of the clock skew 
measurement apparatus and the clock skew measurement method of the present 
invention, a clock skew between clock signals each having a different frequency 

20 from one another can be obtained. By this, a clock signal having relatively low 
frequency can be used as a reference clock signal in a clock skew test, and 
therefore an efficiency and a practical usability of a clock skew test can greatly 
be improved. 

In addition, according to the embodiments of the clock skew 
25 measurement apparatus and the clock skew measurement method of the present 
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invention, a clock skew between clock signals under measurement can be 
estimated by simultaneously sampling in the sequential order clock waveforms 
under measurement utilizing an apparatus that can simultaneously measure two 
channels. By this, the number of required simultaneous samplings can be 
5 decreased from N C2(=N(N-l)/2) to (N- 1) two-channel simultaneous 

measurements, and therefore a measurement time of clock skew can greatly be 
decreased. 

In addition, according to the clock skew measurement apparatus and the 
clock skew measurement method of the present invention, only the minimum 
1 0 number of pins are required when clock signals to be distributed within, for 
example, a semiconductor chip are taken out to the outside of the chip. 
Therefore, the test expense of VLSI test can be reduced. 



